The genome of the blind bee louse fly reveals deep convergences with its social host and illuminates Drosophila origins.

Social insects' nests harbor intruders known as inquilines,1 which are usually related to their hosts.2,3 However, distant non-social inquilines may also show convergences with their hosts,4,5 although the underlying genomic changes remain unclear. We analyzed the genome of the wingless and blind bee louse fly Braula coeca, an inquiline kleptoparasite of the western honey bee, Apis mellifera.6,7 Using large phylogenomic data, we confirmed recent accounts that the bee louse fly is a drosophilid8,9 and showed that it had likely evolved from a sap-breeder ancestor associated with honeydew and scale insects' wax. Unlike many parasites, the bee louse fly genome did not show significant erosion or strict reliance on an endosymbiont, likely due to a relatively recent age of inquilinism. However, we observed a horizontal transfer of a transposon and a striking parallel evolution in a set of gene families between the honey bee and the bee louse fly. Convergences included genes potentially involved in metabolism and immunity and the loss of nearly all bitter-tasting gustatory receptors, in agreement with life in a protective nest and a diet of honey, pollen, and beeswax. Vision and odorant receptor genes also exhibited rapid losses. Only genes whose orthologs in the closely related Drosophila melanogaster respond to honey bee pheromone components or floral aroma were retained, whereas the losses included orthologous receptors responsive to the anti-ovarian honey bee queen pheromones. Hence, deep genomic convergences can underlie major phenotypic transitions during the evolution of inquilinism between non-social parasites and their social hosts.

Dominance mechanisms in supergene alleles controlling butterfly wing pattern variation: insights from gene expression in Heliconius numata.

Loci under balancing selection, where multiple alleles are maintained, offer a relevant opportunity to investigate the role of natural selection in shaping genetic dominance: the high frequency of heterozygotes at these loci has been shown to enable the evolution of dominance among alleles. In the butterfly Heliconius numata, mimetic wing color variations are controlled by an inversion polymorphism of a circa 2 Mb genomic region (supergene P), with strong dominance between sympatric alleles. To test how differences in dominance observed on wing patterns correlate with variations in expression levels throughout the supergene region, we sequenced the complete transcriptome of heterozygotes at the prepupal stage and compared it to corresponding homozygotes. By defining dominance based on non-overlapping ranges of transcript expression between genotypes, we found contrasting patterns of dominance between the supergene and the rest of the genome; the patterns of transcript expression in the heterozygotes were more similar to the expression observed in the dominant homozygotes in the supergene region. Dominance also differed among the three subinversions of the supergene, suggesting possible epistatic interactions among their gene contents underlying dominance evolution. We found the expression pattern of the melanization gene cortex located in the P-region to predict wing pattern phenotype in the heterozygote. We also identify new candidate genes that are potentially involved in mimetic color pattern variations highlighting the relevance of transcriptomic analyses in heterozygotes to pinpoint candidate genes in non-recombining regions.

Genome Size Variation of Chagas Disease Vectors of the Rhodniini Tribe.

The genome size of five Rhodnius species (R. milesi, R. nasutus, R. neivai, R. prolixus, and R. robustus) and two Psammolestes species (P. coroedes and P. tertius) were estimated using flow cytometry and/or k-mer distributions in genome sequences. Phylogenetic generalized linear mixed models highlighted significant genome size variations among species and between sexes, with R. prolixus showing the largest genome. In this study we provide the first data on female genome size in Triatominae. For five species, female genome size did not differ from males, except for R. robustus, where females had smaller genomes. Genome size estimations based on the k-mer distribution method were less than those estimated from flow cytometry, but both methods exhibited the same pattern of sexual differences. Further genomic studies are needed to infer whether genome size variation could be an adaptive trait in Rhodnius.

The fate of a suppressed X-linked meiotic driver: experimental evolution in Drosophila simulans.

Sex-ratio (SR) meiotic drivers are X-linked selfish genetic elements that promote their own transmission by preventing the production of Y-bearing sperm, which usually lowers male fertility. The spread of SR drivers in populations is expected to trigger the evolution of unlinked drive suppressors, a theoretically predicted co-evolution that has been observed in nature. Once completely suppressed, the drivers are expected either to decline if they still affect the fitness of their carriers, or to evolve randomly and possibly get fixed if the suppressors eliminate their deleterious effects. To explore this issue, we used the Paris sex-ratio system of Drosophila simulans in which drive results from the joint effect of two elements on the X chromosome: a segmental duplication and a deficient allele of the HP1D2 gene. We set up six experimental populations starting with 2/3 of X chromosomes carrying both elements (XSR) in a fully suppressing background. We let them evolve independently during almost a hundred generations under strong sexual competition, a condition known to cause the rapid disappearance of unsuppressed Paris XSR in previous experimental populations. In our study, the fate of XSR chromosomes varied among populations, from extinction to their maintenance at a frequency close to the starting one. While the reasons for these variable outcomes are still to be explored, our results show that complete suppression can prevent the demise of an otherwise deleterious XSR chromosome, turning a genetic conflict into cooperation between unlinked loci. Observations in natural populations suggest a contrasting fate of the two elements: disappearance of the duplication and maintenance of deficient HP1D2 alleles.

Evolution of assortative mating following selective introgression of pigmentation genes between two Drosophila species.

Adaptive introgression is ubiquitous in animals, but experimental support for its role in driving speciation remains scarce. In the absence of conscious selection, admixed laboratory strains of Drosophila asymmetrically and progressively lose alleles from one parental species and reproductive isolation against the predominant parent ceases after 10 generations. Here, we selectively introgressed during 1 year light pigmentation genes of D. santomea into the genome of its dark sibling D. yakuba, and vice versa. We found that the pace of phenotypic change differed between the species and the sexes and identified through genome sequencing common as well as distinct introgressed loci in each species. Mating assays showed that assortative mating between introgressed flies and both parental species persisted even after 4 years (~60 generations) from the end of the selection. Those results indicate that selective introgression of as low as 0.5% of the genome can beget morphologically distinct and reproductively isolated strains, two prerequisites for the delimitation of new species. Our findings hence represent a significant step toward understanding the genome-wide dynamics of speciation-through-introgression.

Genome assembly of 3 Amazonian Morpho butterfly species reveals Z-chromosome rearrangements between closely related species living in sympatry.

The genomic processes enabling speciation and species coexistence in sympatry are still largely unknown. Here we describe the whole-genome sequencing and assembly of 3 closely related species from the butterfly genus Morpho: Morpho achilles (Linnaeus, 1758), Morpho helenor (Cramer, 1776), and Morpho deidamia (Höbner, 1819). These large blue butterflies are emblematic species of the Amazonian rainforest. They live in sympatry in a wide range of their geographical distribution and display parallel diversification of dorsal wing color pattern, suggesting local mimicry. By sequencing, assembling, and annotating their genomes, we aim at uncovering prezygotic barriers preventing gene flow between these sympatric species. We found a genome size of  480 Mb for the 3 species and a chromosomal number ranging from 2n = 54 for M. deidamia to 2n = 56 for M. achilles and M. helenor. We also detected inversions on the sex chromosome Z that were differentially fixed between species, suggesting that chromosomal rearrangements may contribute to their reproductive isolation. The annotation of their genomes allowed us to recover in each species at least 12,000 protein-coding genes and to discover duplications of genes potentially involved in prezygotic isolation like genes controlling color discrimination (L-opsin). Altogether, the assembly and the annotation of these 3 new reference genomes open new research avenues into the genomic architecture of speciation and reinforcement in sympatry, establishing Morpho butterflies as a new eco-evolutionary model.

A Population Genomic Assessment of Three Decades of Evolution in a Natural Drosophila Population.

Population genetics seeks to illuminate the forces shaping genetic variation, often based on a single snapshot of genomic variation. However, utilizing multiple sampling times to study changes in allele frequencies can help clarify the relative roles of neutral and non-neutral forces on short time scales. This study compares whole-genome sequence variation of recently collected natural population samples of Drosophila melanogaster against a collection made approximately 35 years prior from the same locality-encompassing roughly 500 generations of evolution. The allele frequency changes between these time points would suggest a relatively small local effective population size on the order of 10,000, significantly smaller than the global effective population size of the species. Some loci display stronger allele frequency changes than would be expected anywhere in the genome under neutrality-most notably the tandem paralogs Cyp6a17 and Cyp6a23, which are impacted by structural variation associated with resistance to pyrethroid insecticides. We find a genome-wide excess of outliers for high genetic differentiation between old and new samples, but a larger number of adaptation targets may have affected SNP-level differentiation versus window differentiation. We also find evidence for strengthening latitudinal allele frequency clines: northern-associated alleles have increased in frequency by an average of nearly 2.5% at SNPs previously identified as clinal outliers, but no such pattern is observed at random SNPs. This project underscores the scientific potential of using multiple sampling time points to investigate how evolution operates in natural populations, by quantifying how genetic variation has changed over ecologically relevant timescales.

Convergent morphology and divergent phenology promote the coexistence of Morpho butterfly species.

The coexistence of closely-related species in sympatry is puzzling because ecological niche proximity imposes strong competition and reproductive interference. A striking example is the widespread wing pattern convergence of several blue-banded Morpho butterfly species with overlapping ranges of distribution. Here we perform a series of field experiments using flying Morpho dummies placed in a natural habitat. We show that similarity in wing colour pattern indeed leads to interspecific territoriality and courtship among sympatric species. In spite of such behavioural interference, demographic inference from genomic data shows that sympatric closely-related Morpho species are genetically isolated. Mark-recapture experiments in the two most closely-related species unravel a strong temporal segregation in patrolling activity of males. Such divergence in phenology reduces the costs of reproductive interference while simultaneously preserving the benefits of convergence in non-reproductive traits in response to common ecological pressures. Henceforth, the evolution of multiple traits may favour species diversification in sympatry by partitioning niche in different dimensions.

Mutation load at a mimicry supergene sheds new light on the evolution of inversion polymorphisms.

Chromosomal inversions are ubiquitous in genomes and often coordinate complex phenotypes, such as the covariation of behavior and morphology in many birds, fishes, insects or mammals1-11. However, why and how inversions become associated with polymorphic traits remains obscure. Here we show that despite a strong selective advantage when they form, inversions accumulate recessive deleterious mutations that generate frequency-dependent selection and promote their maintenance at intermediate frequency. Combining genomics and in vivo fitness analyses in a model butterfly for wing-pattern polymorphism, Heliconius numata, we reveal that three ecologically advantageous inversions have built up a heavy mutational load from the sequential accumulation of deleterious mutations and transposable elements. Inversions associate with sharply reduced viability when homozygous, which prevents them from replacing ancestral chromosome arrangements. Our results suggest that other complex polymorphisms, rather than representing adaptations to competing ecological optima, could evolve because chromosomal rearrangements are intrinsically prone to carrying recessive harmful mutations.

A Variable Genetic Architecture of Melanic Evolution in Drosophila melanogaster.

Unraveling the genetic architecture of adaptive phenotypic divergence is a fundamental quest in evolutionary biology. In Drosophila melanogaster, high-altitude melanism has evolved in separate mountain ranges in sub-Saharan Africa, potentially as an adaptation to UV intensity. We investigated the genetic basis of this melanism in three populations using a new bulk segregant analysis mapping method. We identified 19 distinct QTL regions from nine mapping crosses, with several QTL peaks overlapping between two or all populations, and yet different crosses involving the same melanic population commonly yielded distinct QTL. The strongest QTL often overlapped well-known pigmentation genes, but we typically did not find wide signals of genetic differentiation (FST) between lightly and darkly pigmented populations at these genes. Instead, we found small numbers of highly differentiated SNPs at the probable causative genes. A simulation analysis showed that these patterns of polymorphism were consistent with selection on standing genetic variation. Overall, our results suggest that, even for potentially simpler traits like pigmentation, the complexity of adaptive trait evolution poses important challenges for QTL mapping and population genetic analysis.

The pdm3 Locus Is a Hotspot for Recurrent Evolution of Female-Limited Color Dimorphism in Drosophila.

Sex-limited polymorphisms are an intriguing form of sexual dimorphism that offer unique opportunities to reconstruct the evolutionary changes that decouple male and female traits encoded by a shared genome. We investigated the genetic basis of a Mendelian female-limited color dimorphism (FLCD) that segregates in natural populations of more than 20 species of the Drosophila montium subgroup. In these species, females have alternative abdominal color morphs, light and dark, whereas males have only one color morph in each species. A comprehensive molecular phylogeny of the montium subgroup supports multiple origins of FLCD. Despite this, we mapped FLCD to the same locus in four distantly related species-the transcription factor POU domain motif 3 (pdm3), which acts as a repressor of abdominal pigmentation in D. melanogaster. In D. serrata, FLCD maps to a structural variant in the first intron of pdm3; however, this variant is not found in the three other species-D. kikkawai, D. leontia, and D. burlai-and sequence analysis strongly suggests the pdm3 alleles responsible for FLCD originated independently at least three times. We propose that cis-regulatory changes in pdm3 form sexually dimorphic and monomorphic alleles that segregate within species and are preserved, at least in one species, by structural variation. Surprisingly, pdm3 has not been implicated in the evolution of sex-specific pigmentation outside the montium subgroup, suggesting that the genetic paths to sexual dimorphism may be constrained within a clade but variable across clades.

Recurrent specialization on a toxic fruit in an island Drosophila population.

Recurrent specialization on similar host plants offers a unique opportunity to unravel the evolutionary and genetic mechanisms underlying dietary shifts. Recent studies have focused on ecological races belonging to the same species, but it is hard in many cases to untangle the role of adaptive introgression versus distinct mutations in facilitating recurrent evolution. We discovered on the island of Mayotte a population of the generalist fly Drosophila yakuba that is strictly associated with noni (Morinda citrifolia). This case strongly resembles Drosophila sechellia, a genetically isolated insular relative of D. yakuba whose intensely studied specialization on toxic noni fruits has always been considered a unique event in insect evolution. Experiments revealed that unlike mainland D. yakuba strains, Mayotte flies showed strong olfactory attraction and significant toxin tolerance to noni. Island females strongly discriminated against mainland males, suggesting that dietary adaptation has been accompanied by partial reproductive isolation. Population genomic analysis indicated a recent colonization (∼29 kya), at a time when year-round noni fruits may have presented a predictable resource on the small island, with ongoing migration after colonization. This relatively recent time scale allowed us to search for putatively adaptive loci based on genetic variation. Strong signals of genetic differentiation were found for several detoxification genes, including a major toxin tolerance locus in D. sechellia Our results suggest that recurrent evolution on a toxic resource can involve similar historical events and common genetic bases, and they establish an important genetic system for the study of early stages of ecological specialization and speciation.

Ancient balancing selection at tan underlies female colour dimorphism in Drosophila erecta.

Dimorphic traits are ubiquitous in nature, but the evolutionary factors leading to dimorphism are largely unclear. We investigate a potential case of sexual mimicry in Drosophila erecta, in which females show contrasting resemblance to males. We map the genetic basis of this sex-limited colour dimorphism to a region containing the gene tan. We find a striking signal of ancient balancing selection at the 'male-specific enhancer' of tan, with exceptionally high sequence divergence between light and dark alleles, suggesting that this dimorphism has been adaptively maintained for millions of years. Using transgenic reporter assays, we confirm that these enhancer alleles encode expression differences that are predicted to generate this pigmentation dimorphism. These results are compatible with the theoretical prediction that divergent phenotypes maintained by selection can evolve simple genetic architectures.

Pigmentation in Drosophila melanogaster reaches its maximum in Ethiopia and correlates most strongly with ultra-violet radiation in sub-Saharan Africa.

BACKGROUND: Pigmentation has a long history of investigation in evolutionary biology. In Drosophila melanogaster, latitudinal and altitudinal clines have been found but their underlying causes remain unclear. Moreover, most studies were conducted on cosmopolitan populations which have a relatively low level of genetic structure and diversity compared to sub-Saharan African populations. We investigated: 1) the correlation between pigmentation traits within and between the thorax and the fourth abdominal segment, and 2) their associations with different geographical and ecological variables, using 710 lines belonging to 30 sub-Saharan and cosmopolitan populations. RESULTS: Pigmentation clines substantially differed between sub-Saharan and cosmopolitan populations. While positive correlations with latitude have previously been described in Europe, India and Australia, in agreement with Bogert's rule or the thermal melanism hypothesis, we found a significant negative correlation in Africa. This correlation persisted even after correction for altitude, which in its turn showed a positive correlation with pigmentation independently from latitude. More importantly, we found that thoracic pigmentation reaches its maximal values in this species in high-altitude populations of Ethiopia (1,600-3,100 m). Ethiopian flies have a diffuse wide thoracic trident making the mesonotum and the head almost black, a phenotype that is absent from all other sub-Saharan or cosmopolitan populations including high-altitude flies from Peru (~3,400 m). Ecological analyses indicated that the variable most predictive of pigmentation in Africa, especially for the thorax, was ultra-violet (UV) intensity, consistent with the so-called Gloger's rule invoking a role of melanin in UV protection. CONCLUSION: Our data suggest that different environmental factors may shape clinal variation in tropical and temperate regions, and may lead to the evolution of different degrees of melanism in different high altitude populations in the tropics.

Local dynamics of a fast-evolving sex-ratio system in Drosophila simulans.

By distorting Mendelian transmission to their own advantage, X-linked meiotic drive elements can rapidly spread in natural populations, generating a sex-ratio bias. One expected consequence is the triggering of a co-evolutionary arms race between the sex chromosome that carries the distorter and suppressors counteracting its effect. Such an arms race has been theoretically and experimentally established and can have many evolutionary consequences. However, its dynamics in contemporary populations is still poorly documented. Here, we investigate the fate of the young X-linked Paris driver in Drosophila simulans from sub-Saharan Africa to the Middle East. We provide the first example of the early dynamics of distorters and suppressors: we find consistent evidence that the driving chromosomes have been rising in the Middle East during the last decade. In addition, identical haplotypes are at high frequencies around the two co-evolving drive loci in remote populations, implying that the driving X chromosomes share a recent common ancestor and suggesting that East Africa could be the cradle of the Paris driver. The segmental duplication associated with drive presents an unusual structure in West Africa, which could reflect a secondary state of the driver. Together with our previous demonstration of driver decline in the Indian Ocean where suppression is complete, these data provide a unique picture of the complex dynamics of a co-evolutionary arms race currently taking place in natural populations of D. simulans.

A genome-wide, fine-scale map of natural pigmentation variation in Drosophila melanogaster.

Various approaches can be applied to uncover the genetic basis of natural phenotypic variation, each with their specific strengths and limitations. Here, we use a replicated genome-wide association approach (Pool-GWAS) to fine-scale map genomic regions contributing to natural variation in female abdominal pigmentation in Drosophila melanogaster, a trait that is highly variable in natural populations and highly heritable in the laboratory. We examined abdominal pigmentation phenotypes in approximately 8000 female European D. melanogaster, isolating 1000 individuals with extreme phenotypes. We then used whole-genome Illumina sequencing to identify single nucleotide polymorphisms (SNPs) segregating in our sample, and tested these for associations with pigmentation by contrasting allele frequencies between replicate pools of light and dark individuals. We identify two small regions near the pigmentation genes tan and bric-à-brac 1, both corresponding to known cis-regulatory regions, which contain SNPs showing significant associations with pigmentation variation. While the Pool-GWAS approach suffers some limitations, its cost advantage facilitates replication and it can be applied to any non-model system with an available reference genome.

Rapid rise and fall of selfish sex-ratio X chromosomes in Drosophila simulans: spatiotemporal analysis of phenotypic and molecular data.

Sex-ratio drive, which has been documented in several Drosophila species, is induced by X-linked segregation distorters. Contrary to Mendel's law of independent assortment, the sex-ratio chromosome (X(SR)) is inherited by more than half the offspring of carrier males, resulting in a female-biased sex ratio. This segregation advantage allows X(SR) to spread in populations, even if it is not beneficial for the carriers. In the cosmopolitan species D. simulans, the Paris sex-ratio is caused by recently emerged selfish X(SR) chromosomes. These chromosomes have triggered an intragenomic conflict, and their propagation has been halted over a large area by the evolution of complete drive suppression. Previous molecular population genetics analyses revealed a selective sweep indicating that the invasion of X(SR) chromosomes was very recent in Madagascar (likely less than 100 years ago). Here, we show that X(SR) chromosomes are now declining at this location as well as in Mayotte and Kenya. Drive suppression is complete in the three populations, which display little genetic differentiation and share swept haplotypes, attesting to a common and very recent ancestry of the X(SR) chromosomes. Patterns of DNA sequence variation also indicate a fitness cost of the segmental duplication involved in drive. The data suggest that X(SR) chromosomes started declining first on the African continent, then in Mayotte, and finally in Madagascar and strongly support a scenario of rapid cycling of X chromosomes. Once drive suppression has evolved, standard X(ST) chromosomes locally replace costly X(SR) chromosomes in a few decades.

Molecular signature of epistatic selection: interrogating genetic interactions in the sex-ratio meiotic drive of Drosophila simulans.

Fine scale analyses of signatures of selection allow assessing quantitative aspects of a species' evolutionary genetic history, such as the strength of selection on genes. When several selected loci lie in the same genomic region, their epistatic interactions may also be investigated. Here, we study how the neutral polymorphism pattern was shaped by two close recombining loci that cause 'sex-ratio' meiotic drive in Drosophila simulans, as an example of strong selection with potentially strong epistasis. We compare the polymorphism data observed in a natural population with the results of forward stochastic simulations under several contexts of epistasis between the candidate loci for the drive. We compute the likelihood of different possible scenarios, in order to determine which configuration is most consistent with the data. Our results highlight that fine scale analyses of well-chosen candidate genomic regions provide information-rich data that can be used to investigate the genotype-phenotype-fitness map, which can hardly be studied in genome-wide analyses. We also emphasize that initial conditions and time of observation (here, time after the interruption of a partial selective sweep) are crucial parameters in the interpretation of real data, while these are often overlooked in theoretical studies.